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M.N.  Woinsky  proposed  the  application  of  various  non-parametric 
statistics  to  spectral  data  in  order  to  detect  narrow  band  signals. 

He  determined  the  asymptotic  relative  efficiency  (ARE)  of  two  such  tests. 

Hansen  and  Olsen  have  found  the  distribution  of  ranks  of  spectral  data 
for  two  types  of  signal  embedded  in  noise.  Here  the  distribution  of  ranks 
for  a third  type  of  signal  is  derived.  The  performance  of  various  non- 

parametric  statistics  for  these  three  signals  are  calculated.  j 
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I.  INTRODUCTION 

Following  Woinsky  \^1^  we  consider  the  situation  in  which 
a signal  embedded  in  a time  series  appears  as  a narrow  line  in  the 
frequency  domain.  Detection  is  based  on  the  matrix  of  spectral  estimates 
j : i*l,2,...,M;  j=l,2,...,N.  Here  time  is  referenced  by  the  i parameter 
and  frequency  by  the  j parameter.  The  parameter  N is  the  largest  number 
of  frequency  estimates  for  which  the  spectrum  can  be  said  to  be  flat,  while 
M is  the  number  of  spectral  estimates  per  frequency  cell  to  be  used  in 


the  decision  process. 

Let  R be  the  rank  of  the  for  a fixed  time  interval, 
that  is  i fixed,  j*l,2,...,N.  The  null  hypothesis  is  that  all  N frequencies 
contain  noise  only.  For  the  X independent  identically  distributed  random 
variables,  the  R are  uniformly  distributed  over  the  integers  1,2,...N.  The 
alternative  hypothesis  is  that  a signal  occurs  at  one  of  the  frequencies 


represented.  In  order  t?o  determine  the  distribution  of  the  ranks  in 


this  case,  one  must  assume  a specific  distribution  for  the  ; that  is, 
one  must  specify  the  nature  of  the  signal  and  the  noise. 


II.  DISTRIBUTION  OF  THE  RANK  STATISTICS 
Suppose  X^.X^,...,^  represent  a set  of  N independent  random 
variable  spectral  estimates  of  which  the  first  N-l  correspond  to  frequencies 

V 

without  signal  components  and  of  which  the  last,  XN»  corresponds  to  a 
frequency  cell  at  which  a signal  occurs.  is  denoted  as  the  signal  random 
variable.  We  consider  the  case  where  the  first  N-l  estimates  are  exponen- 
tial random  variables  with  mean  that  is,  the  time-bandwidth  product 
for  the  assumed  spectral  analyzer  is  unity  and  ¥ represents  the  expected 
power  in  each  cell.  This  is  a reasonable  assumption  as  theoretical  , 12] 
and  empirical  studies  ^3^  have  shown  noise  spectral  estimates  to  be  independent 
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I 


J 


exponential  random  variables  for  a variety  of  time  series. 
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Three  types  of  signal  are  considered  below.  Case  I refers  to 
an  X^  which  is  an  exponential  random  variable  with  parameter  0(1+V);  that  is, 
the  signal  is  a "Gaussian  target"  and  0 is  cell  signal-to-noise  power  ratio  CO- 
Case  II  refers  to  a distribution  of  3^  defined  by 

f(V  = 1 exp(“  K1  XN46)  *o  [<2*'lexN)!5]  ••••<!) 


where  x„  > 0 and  I (x)  is  the  modified  Bessel  function  of  order  zero. 

N — o 

This  distribution  occurs  when  the  signal  acts  like  a sinusoidal  target  10  0]  • 
Case  III,  formulated  by  Swerling  for  a model  of  a fluctuating 
radar  target  OHO  , refers  to  a density  function  of  the  form 

f(xN)  = \ bV1(l  + axNy_1)  exp(-  \ x^f'1)  (2) 


where 


b represents  (1  + -^6)  * 


represents  L2  + (4/0)"]  _1.  1 


Let  P(k)  be  the  probability  that  the  signal  frequency  takes  rank  k, 
where  k = 1,2,...N.  This  occurs  if  there  are  k-1  noise  random  variables 
less  than  and  N-k  such  random  variables  greater  than  X^. 

We  find  the  probability  of  one  such  ordering  and  then  multiply 


J. 

by  for  the  final  expression  for  P(k) . 

the  probability  of  one  specific  ordering  is 


since  are  independent. 


[k-1  specified  noise  fre- 
quency spectral  estimates 
are  less  than  x 


remaining  N-k 
estimates  are  dx...(3) 
greater  than  x 


where  P denotes  the  probability  of  an  event  and  f(x)  is  the  density  function 


(a)  CASE  I 

In  Case  I,  expression  (3)  Cakes  the  form 


f 2*1' u+0)  expr  ^YCl+G) 

o 


e , 

f -i  . x k-l  « I N-k 

l"  /2y(i+©)J  — [ exp(“  y/2v)dy  I “ [ exp(-  y/2y)dy  I dx 


2V  (1+0)  fQ  1-exP("  ! 2 ^ 


! k-i  r 


exp(-  / ^ exp  (-  / 


2V  (1+0) ' 


-at  I |_l-exp(-y)J  k 1 [exp(-y)]  N k exp(-oy)dy, 

o 

where  a represents  (1+6)  1 and  y represents  x/(2V). 
Now,  by  [8]  page  305  Formula  3.312-1, 


00  -J 

/ | 1-exp (-  X/g)|  V 1 exp(-px)dx  = 8B(8u,vy 

if  Re  8>o,  Re  v>o.  Re  p>o,  where  B (x,y)  represents  the  Beta  function.* 

Thus  our  integral  reduces  to 

aB(N-k  + a ,k)  - aT(N-k+a)r(k)/r(N+a) 

We  note  that  the  conditions  on  8 and  v are  always  satisfied;  care  must  be 
taken  with  the  condition  on  p which  reduces  to  the  requirement  that 
N+a>k.  This  holds  if  a>0.  If  a«0  we  have  the  case  where  cell  signal-to-noise 
power  ratio  is  infinite.  It  can  be  independently  shown  that  the  same 
result  holds  in  this  case. 

Finally,  multiplication  by  (^_^)  gives  the  result 


P(k)  - P(k  N a)  - arW.r(N+u-k)  ...(4) 
rite, N, a;  r (N+a)  T (N-k+1) 

for  k = 1....N.  This  equation  is  equivalent,  although  simpler  in  form,  to 
equation  (32)  of  G). 

It  is  easily  seen  that  if  0 is  zero,  P(k)  is  N * for  all  k. 


the  uniform  distribution.  The  ranks  are  uniformly  distributed  for  any 
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distribution  of  the  as  long  as  the  X^  are  independent  and  identically 
distributed . 

It  can  be  shown  that  P(k)  <_  P(k')  for  k < k',  with  equality 
holding  only  if  6 equals  zero. 


(b)  CASE  II 

In  this  case  expression  (3)  becomes  the  integral 

!J-k 


/0  [l-exp(-y)]  k 1 j^exp(-y)  jN  k exp  -(y+0)  ] Iq  | 2(y6)'5j  dy 

= exp(-0)  /o  £l-exp(-y)|  k 1 exp  | -y(N-k+l)  | IQ^2(y0)^ 

k“  X 00  ^ «•  *■ 

- exp(-0)£  (-l)1(k,1)  / exp  f-y(N-k+i+l)]l  1 2(y0)'5 

i=o  i o t J ° 


dy 

dy 


where  y represents  x 4'-1. 


Now,  by  [j.0^]  Pa8e  Formulae  29.3.81, 

oo  r i ■y 

jQ  exp(-pt)IQ  I 2(at)^ J dt  = p exp(a/p),  where  Re  p>0. 


Hence  our  integral  becomes 
k-1 


exp(-0)  l (k“X)  (-1)1  N-_-kTi4-T  exp  [ 0/(N-k+i+l] 
i=o  ' J 

whenever  Re  N-k+i+k>0.  This  condition  is  always  satisfied. 


,N-1, 


Multiplying  by  (^_^)  yields,  after  simplifying, 

r(N)exp(-0)  7 exP.  [oj  (N+l-k+i)|  . . 

l.  (N+i-k+i)ra+nra-i) 


P(k)  = P ( k , N , 0) 


T(N-k+l) 


i=o 
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for  k “ 1,2, . . .N. 

This  is  identical  to  equation  (30)  of  • 

Equation  (5)  gives  the  density  function  for  Case  II.  Care  must 
be  taken  in  its  evaluation,  as  we  are  summing  a strictly  alternating  series 
whose  sum  while  lying  between  zero  and  one  contains  terms  of  different  orders 
of  magnitude.  It  can  be  shown  analytically  that 

P(k)  ■ -s  if  6 > 0. 


(c)  CASE  III 

For  the  density  function  given  by  Equation  (2),  expression 
(3)  becomes 

\ Jo  bW+ax*"1)  exp(-  -Ixb'T1)  [ exp(-  N_k 

| 1-exp (-  -IxH1  "*■)  Jk  ^dx 


. 2 Vr  f , , 1 

k-1  f "1 

- b ) / l-exp(-y) 

exp  -y(N-k+b) 

> o ^ 

L * 

00  r 

+ 2a  JQ  £ l-exp(-y)j 

k 1 exp  £ -y(N-k+b)j 

where  y represents  -IxH1  ^ . 


oo 

J0  [_l-exp(-y)J  k 1 exp  ^-y(N-k+b)J  dy 


- B(N+k+b,k) 

* r(N-k+b)r(k)/r(N+b). 

Also,  by  M formula  5.9, 

/o  ^ l-exp(-y)  ] k 1 exp  £-y  (N-k+b)  ] ydy 
= B(N-k+b,k)  £ D(N+b)  - D(N-k+b)*| 

where  D represents  the  Digamma  function.  From  [_io]  formula  6.36  the 

k~i 

expression  in  the  square  bracket  equals  £ (m+N-k+b)  , multiplication  by 

n-1  m“° 

(J.J)  Yields 


P(k)  = P(k,N,b)  = b2  (N-k+b)  Q( N+b ) (N-k+l)l  _1 

T k_1  il 

1 + 2a  l (m+N+b-k) 
m=o 

As  inCase  I it  can  be  shown  analytically  that  P(k)  £ P(k') 
for  k < k'  with  equality  holding  only  if  6 equals  zero. 

Figure  1 contains  typical  density  functions  when  N = 8.  For 
low  cell  signal-to-noise  ratios  these  functions  are  almost  the  same 
(Fig. 1(a)),  whereas  for  larger  ratios  that  of  Case  II  tends  to  the  limiting 
case  P(N)  = 1 "faster"  than  the  density  function^jf  Case  I.  The  results 
for  Case  III  are  intermediate,  and  are  not  included  in  the  Figures. 


III.  DETECTOR  PERFORMANCE 

The  modified  nonparametric  decision  scheme  is  based  fin  the  sum  of 
M independent  random  variables  representing  some  function  of  the  rank  at  each 
frequency.  That  is,  we  decide  a signal  is  present  at  frequency  j depending 


on  the  value  of 


5j  " M ^=ia(Rij)  •" 


where  a(x)  is  some  function  of  x.  For  M small  we  use  numerical  convolution, 

for  M large  we  apply  the  Central  Limit  Theorem  to  obtain  the  density  function  of  S. 

We  determine  S . when:  * 

J 

(1)  a(R  ) = R (modified  Mann-Whitney  statistic). 


ij'  “ij 


(2)  a(R  ) 


l (modified  Savage  test). 


£=N-R . .+1 
ij 


(3)  a(R„)  = sign(Rij  - N/2)  (modified  Median  test). 

Dillard  to  has  studied  case  a(Rij)=  sign  (R^-N) . The  results 
above  would  enable  us  to  study  the  detection  a(Rij)  = sign  (R^-kN)  for 
0 < k < 1.  In  these  cases  of  course  tables  of  the  binomial  distribution 


may  be  used  to  determine  false  alarm  and  detection  probabilities 


K 
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IV  RESULTS 

Computer  programs  were  written  to  determine  the  performance 
under  various  combinations  of  N,  M and  0 of  the  three  detectors  above  in 
Cases  I,  II  and  III. 

Figure  2 gives  some  typical  results.  In  most  cases  we  present 
the  Savage  statistic  results,  the  "best"  detector  considered  here. 

In  all  cases  considered,  the  Savage  detector  gave  better 
detection  results  than  the  Mann-Whitney  detector.  In  turn  the  Mann- 
Whitney  performed  better  than  the  Median  detector.  Woinsky  • L12) 

has  shown  with  regard  to  the  ARE  that  the  Savage  statistic  performs  better 
than  the  Mann-Whitney . Figure  2(a)  gives  a typical  result  for  a Case  II 
situation  with  low  signal-to-noise  ratio  and  a large  sample  size.  Case  I, 
II  and  III  signals  have  the  same  performance  for  the  various  detectors  at 
very  low  signal-to-noise  ratios,  at  other  values  for  a fixed  false  alarm 
rate  we  can  detect  a sinusoidal  target  with  greatest  probability.  Figures 
2(b)  and  (c)  illustrate  this. 

V.  CONCLUSIONS 

We  have  determined  the  density  functions  of  the  rank  statistic 
for  three  realistic  signal  models,  hence  exact  receiver  performance  can 
be  evaluated  for  any  practical  detector  based  on  these  ranks. 

Illustrative  numerical  results  were  presented  for  three  weltl- 

I 

known  rank  statistics. 

The  result  of  Woinsky  indicating  the  superiority  of  the  Savage 
statistic  over  the  Mann-Whitney  statistic  in  the  ARE  sense  has  been  found 
applicable  to  all  cases  tested  through  numerical  evaluation  of  the  exact 
expressions  derived  over  a wide  range  of  signal-to-noise  ratio  and  sample 
size. 


# 
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ACCUMULANTS 


Figure  11(c).  Accumulants  versus  probability  of  detection  for  sinusoidal 
and  Gaussian  signals.  Savage  Test;  16  frequencies  ranked; 
cell  signal-to-noise  ratio  -6  dB,  false  alarm  rate  10“^. 
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